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Abstract: It is well known that the renormalization group equations depend on the scale 
where they are applied. This phenomenon is especially relevant for the massive fields in 
curved space, because the decoupling effects may be responsible for important cosmological 
applications like the graceful exit from the inflation and low-energy quantum dynamics of 
the cosmological constant. We investigate, using both covariant and non-covariant methods 
of calculations and mass-dependent renormalization scheme, the vacuum quantum effects of 
a massive scalar field in curved space-time. In the higher derivative sector we arrive at the 
explicit form of decoupling and obtain the /3-functions in both UV and IR regimes as the 
limits of general expressions. For the cosmological and Newton constants the corresponding 
/^-functions are not accessible in the perturbative regime and in particular the form of 
decoupling remains unclear. 
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1. Introduction 

The renormalization group methods represent one of the most universal tools in the Quan- 
tum Field Theory and beyond: from string theory to lattice QCD and statistical mechanics. 
In the last decades, there was a growing interest to the renormalization group for quan- 
tum fields in curved space-time, especially in relation to the anomalies and cosmological 
applications. The renormalization group and related phenomena of the decoupling of the 
massive quantum matter fields in an external gravitational field are especially relevant for 
such phenomena as the graceful exit from the anomaly-induced inflation [EL ||] and the 
possible low-energy quantum dynamics of the cosmological constant Q. The relevance of 
these applications (which may help in better understanding of the coincidence puzzle for 
the ACDM cosmological model and also may lead to a natural theory of inflation) requires 
a detailed quantitative analysis of the physical aspects of the renormalization group. In 
this paper we are going to start this investigation. 

The standard approach to the renormalization group in an external gravitational field 
has been formulated in [Q, [5| (see also B, |7]]) on the basis of the Minimal Subtraction 
(MS) renormalization scheme or the modified Minimal Subtraction scheme MS. In the 
framework of the MS scheme, the renormalization of quantum matter in curved space- 
time is well understood The principal results are quite simple and look as follows. 
Starting from the renormalizable QFT in the flat space-time, one can always construct the 
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renormalizable QFT in curved space-time. In particular, if we treat the gravitational field 
as a perturbation over the flat background, then the presence of an external gravity does 
not increase the superficial degree of divergence for a given Feynman diagram. Taking the 
covariance and locality considerations into account, one can easily determine the complete 
set of possible divergent structures and identify the corresponding necessary counterterms. 
In general, the action of the renormalizable theory includes the covariant generalization of 
the action in flat space plus a set of curvature-dependent local terms, which are necessary 
for the renormalizability. In particular, one has to introduce a special action of vacuum 
(external gravity 1 ) which has the form 

Svacuuin — Seh + Shd ■ (1.1) 
Here the first term is the Einstein-Hilbert action with cosmological constant. 2 

Seh = ~ J d A xg l l 2 (R + 2A) . (1.2) 

The second action contains necessary higher derivative terms 

Shd = J d A x g l/2 { aiC 2 + a 2 R 2 + a 3 E + a A V 2 R } , (1.3) 



where 

''fj.uap 



C 2 = R 2 p -2R 2 +l/3R 2 (1.4) 



is the square of the Weyl tensor and 

E = RauafsR^ ~ 4 R af sR aP + R 2 (1.5) 

is the integrand of the Gauss-Bonnet topological invariant. ai,...,04 (and also G and A 
in ( |1.2| )) are the parameters of the vacuum action. In the present paper we shall focus our 
attention on the renormalization of the vacuum parameters Gtj, G, A and on the quantum 
corrections to the vacuum action (1.1). 



For the interacting fields there is also renormalization in the matter fields sector. Let 
us remark that the renormalization of those terms that have direct analogs in flat space is 
exactly the same as in the flat space ||]. In addition, there is also the renormalization of 
the nonminimal parameter £ in the scalar action 

Ss = \ f d"xg l l 2 {r^^ + mV^^ 2 } . (1.6) 

The renormalization of £ has many interesting features (see, e.g. |9], |l(|) but is beyond 
the scope of the present paper because, in the cosmological framework, the effect of this 
renormalization is less important than the renormalization of vacuum terms. 



x Let us emphasize that gravity itself is not quantized in this paper, so we are confined to the so-called 
semiclassical approach. 

2 Our notations are Euclidean metric r\^ v = diag(+ + H — h) and the definition of the curvature tensor 

- iiv ■ ■ ■ 
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The vacuum divergences of the theory fll.q ) can be easily calculated using the Schwinger- 
DeWitt method and dimensional regularization 



+5 (« " + " I) * + 5™*} ■ (L7) 

where uj is the parameter of the dimensional regularization (2co is the space-time dimen- 
sion) and fx is a dimensional parameter which is introduced in order to compensate the 
change of the dimension. The renormalization of the vacuum action in the MS scheme 
proceeds as follows. One has to introduce the local counterterm AS and then compare 
the bare So[A(°), G^°\ a\ ] and renormalized Sij[A, G, en] = S + AS actions. This leads to 
the standard renormalization relations (see, e.g. ||). In particular, 



A(°) 
8ttG(°) 



A m 
+ 



i 



8vrG 4(4vr) 2 (u; - 2) 



(1.8) 



Taking into account the fact that A^ ) does not depend on fi, we arrive at the usual (see, 
e.g. [|ll]]) renormalization group equation for the cosmological constant 

where we have set u = 2. Using similar calculations, we arrive at the renormalization 
group equations for other vacuum parameters (see, e.g. Q) 



dfi ^ v ; 120 (4tt) 2 ' 

^ = fi ,<nB> = -^(t-\y. (i.io) 

Let us emphasize that these equations were obtained within covariant formalism and with- 
out imposing restrictions on the external metric g^ v . 

The interpretation of the renormalization group equations in curved space-time re- 
quires a special attention. Of course, this is a usual situation for the MS scheme, where 
even in flat space one has to make a special effort in order to interpret the parameter 
fi. Since /i is an auxiliary parameter introduced in the dimensional regularization, to get 
physical results, one should trade it for some physical quantity. This procedure may lead to 
ambiguities, which can be resolved in the framework of the more physical mass-dependent 
scheme of renormalization. For the usual flat-space theories, say in QED or Yang-Mills, 
the interpretation of fi can be achieved as follows. One can compare the renormalization 
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group in the MS scheme and in a mass-dependent scheme, which must coincide in the high 
energy limit. In this way, \i acquires a natural interpretation. Depending on the situation 
this may be the energy of scattering of quanta, VEV of the scalar field or temperature etc. 

On the other hand, at low energies, the physically correct result requires the use of a 
mass-dependent scheme. But, in the case of the vacuum renormalization in a curved space- 
time, the notion of energy for the gravitational quanta is not well defined, and the relation 
between MS and mass-dependent scheme is a priori unclear. At the same time, there 
is a considerable interest in implementing ideas of effective Quantum Field Theory into 
the gravitational framework: both for the classical gravitational background and effective 
quantum gravity [O]. To this end, we need to achieve a quantitative description of how 
the decoupling of the heavy particles in the external gravitational field occurs. 

Let us briefly remind how decoupling proceeds in flat space in a mass-dependent scheme 
by using the QED example (see, e.g. Hlg] ). The 1-loop vacuum polarization in the dimen- 
sional regularization is equal to 



-.2 



<~- . > r 1 7 f , M s m 2 +px(l — x) 



-^{PvPv-P 9ixv) — — ^-«~/ dxx(l - x) In — - — — ^ , (1.11 



6(w-2) 6 Jo 4vrAt 2 



where p is an external momentum, m is fermion mass, 7 is Euler's constant, and /j, is 
the dimensional parameter of the regularization. By subtracting in the MS scheme the 
^2 -pole an d the 7-term, one find the following contribution: 

e 2 , 2 x f 1 , x , m 2 e + p 2 x{\ - x) , 
- 7^2 (PvPv ~ P 9ttv) ax x(l - x) In ^—^ . (1.12) 

The /3-function in the MS scheme is obtained by calculating the derivative (e/2)nd/dfi 
acting on the coefficient of the operator {PfiPv — P 2 9fj.v)- The result is well known 

ftP) = (1.13) 



Thus, the /3-function in the MS scheme is a constant independent of the fermion mass and 
/x. To calculate the /3-function in a mass-dependent scheme such as, e.g., the momentum 
space subtraction scheme, one has to subtract the value of the graph at a point p 2 = M 2 
and calculate the derivative (e/2)Md/dM acting on the coefficient of the (p^p u —p 2 9^v) 
projector. As a result, we obtain 

e 3 f 1 . M 2 x(l -x) 

d e = ;r / dxx(l-x) 7: — - ; — ; r. (1T4) 

y 2ir 2 J y 1 m 2 + M 2 x(l - x) y ' 

It is easy to verify that when the renormalization point is much larger than the fermion 
mass, the /3-function obtained in the mass-dependent scheme coincides with the MS (3- 
function ( 1.13j ). As the renormalization point passes through m e , massive fermion decouples 



and, for M <C m e , its contribution to the /3-function vanishes as 



e 3 M 2 

/ 9 e (M< A i e « — ^ — . (1.15) 
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Thus, the MS /3-function gives wrong result at energies less than m e . Decoupling of heavy 
particles is implemented in the MS scheme by putting, artificially, the /3-function to zero 
for fj, < m. This can be called "hard decoupling" or "sharp cut-off" approximation. Indeed, 
this is just an approximation to the real picture of decoupling, which may be essentially 
scheme-dependent. In particular, the example of QED shows that, in principle, the decou- 
pling in gravity can proceed in a soft manner such that the /3-function is continuous at the 
scale m. 

The purpose of this paper is to perform, for the first time, the derivation of the 
renormalization group equations for the parameters A, G, a\ and 02 in a simple mass- 
dependent scheme. We consider a massive non-minimally coupled scalar theory, but the 
results can be also generalized for massive fermions and vectors. In section 2 we present 
the derivation of the polarization operator of gravitons using dimensional regularization. 
The table of cumbersome integrals are exposed in the Appendix A. Although we will later 
confirm our principal conclusions by using covariant heat kernel calculations, we believe 
that the linearized gravity calculations are useful. These calculations are performed on the 
flat space background and, hence, they are maximally close to usual perturbative quantum 
field theory calculations, where we have a good intuition of how the decoupling of massive 
particles works. In section 3, the polarization operator of section 2 is verified using the 



previously known expressions for the summation of the Schwinger-DeWitt series [14, 15]. 
In section 4, the renormalization group equations for the vacuum parameters are derived 
and the comparison with similar equations in the MS scheme is given. Another calcula- 
tion which is aimed to clarify the relation between the MS and mass-dependent scheme is 
presented in the Appendix B, where we apply the covariant cut-off regularization of the 
proper-time integral representation for the effective action. 

Unfortunately, in this regularization scheme one can recover the reliable form of the 
(3- functions only at high energy limit, while the small cut-off limit is inconsistent. In section 
5 we discuss the results of the calculations in the mass-dependent scheme. In particular, 
this opportunity is used to discuss the form of the effective action which could correspond 
to the running of the cosmological and inverse Newton constants at low energies. After all, 
we draw our conclusions. 

In order to avoid the confusion, let us fix the notations for the distinct types of the 
/3-functions. We shall denote, in this paper, the MS-scheme /3-function for the effective 
charge C as /3c(MS) and the /3-function for the same object, derived in a mass-dependent 

satisfy the relation 
2. The non-covariant perturbative calculation 

The renormalization group for the parameters of the vacuum action can be established 
through the perturbative calculation of the quantum corrections to the gravitational 2- 



scheme as /3c. The high energy limit of this /3-function will be denoted as /3q V and the 
low-energy limit as (3l R . Of course, we expect that the correctly defined /3-function would 
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point function. Therefore, in this section we shall derive the contributions of the loop 
of massive particle to the propagator of the gravitational perturbations on the flat 
background 



(2.1) 



The calculations will be performed in the Euclidean space. The diagrams which contribute 
to the gravitational polarization operators are presented at Figure 1. 




(b) 





Figure 1. The 1-loop diagrams with the loop of massive scalar and two external 
gravitational lines. The straight lines correspond to the massive scalar and wavy lines 
to the external field h^ v . 



The 3-point and 4-point vertices (see Figure 2) have the form 

V al3 (k,p, q) = (2vr) 4 [^(T] a(3 p ■ q - q a p^ - p a q P ) + t{k a k P - k 2 rf p ) - ^V^] i (2-2) 



V^' a/3 (k,p,q,l) = (2tt) 4 { 5^' al3 [z(k 2 + l 2 + 3k-l)-^p-q 



+ 



[ \p-q-\{k 2 + l 2 + k-l) 



+ rj 



+ 



a/3 



/JO 



p u q< 3 + q v pP - 2i (k v kP + l v f) - 2i k v l? - £1"^ } , (2.3) 



where, for the sake of compactness, we do not maintain the symmetry inside the couples 
of indices (a/3) and (fj-v). 



WW\AA* 

k, a/3 \ 

\ V 



I, afi 



k, \iv 




P 



Figure 2. Two relevant vertices of the gravity-scalar interactions. 
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The first observation is that the diagram (lc) does not contribute to the renormal- 
ization group because we are dealing with the semiclassical approximation. The object 
of interest is the effective action of the external metric, and its expansion into series in 
hn V does not lead to the diagram (lc). For this reason we have to consider only (la) 
and (16) diagrams. Let us start from (la). The lengthy calculations in the dimensional 
regularization give the following expression for the polarization operator 

lC' a/? (p) = "^V^ (p 4 /i-2p 4 / 2 + 6p 2 /3 + 2p 4 I 4 + 8/ 9 ) - ipVWV7-2J 5 +i4)- 

- \ {j) m p v pP + t P P u p a + rp"/ + (4/a - 4Ze + h) - 5^ ' af3 I 9 - 

8 

- I ( VTP V + V al3 P V ) (P 2 h - P 2 h + 4I 8 - 4J 6 ) - f h W ~ p\n (p V - p 2 V al3 ) - 

o Z 

- I (?V - pV") V aP (p 2 h + 4/ 3 ) - I (pV - pV) *T (p 2 A + 4/ 3 ) + 
o o 

+ 7 (// " pV)j>V(/i - 2/4) + 7 (pV " pV' 3 ) pV(Ii " 2I 4 ) + (2.4) 



+ 4 (4vr)- 1 j 71 V ' 



where the integrals I\, I9 are defined in the Appendix A. 

The second diagram is simpler and the corresponding expression is 



n 



16 



\p) 



4 (4vr) 
1 



P 



2 , a/3 



+ 



+ 



r(-a;) (u; - 2) 
8(47r)™ 



(m 2 ) w (2o^' Q/3 -rfrf 1 ^ 



(2.5) 



In order to apply the mass-dependent renormalization scheme, we need both divergent 
and finite parts of the polarization operator. Let us start from divergencies and establish 
the correspondence with the covariant expression (|l.7|). 

For this end, since we are working with the polarization operator, the terms in the 
vacuum action ( |1.1[ ) must be expanded up to the second order in h^ u . As far as the Gauss- 
Bonnet term does not contribute to the propagator, we can replace the C^^-term by the 
expression 2W, where W = R 2 ^ — \R 2 - Then the relevant bilinear expressions have the 
form 

jd^xg 1 ' 2 = Jd^xh^^v^Vap- \^u, a p)h af3 +..., (2.6) 



J d 4 xg^ 2 R = J J 



x hT 



-7- 



h a(i + 



(2.7) 



d 4 xg 1/2 R 2 = I d A x h? v 



h al3 + .. 



(21 



d 4 xq 1/2 W = / d A x h" v 



+ ^ (V^d a df3 + r} a pd^d v ) d 2 



h af3 + .. 



where the dots stand for the lower and higher order terms. 
Taking into account the formulas 



r(2 - u) 
r(i-u) = 



1 



2-uj 
1 



2-cu 



h 0(2-u), 
1 + 0(2-uj), 



(2. 



r(-w) 



' + ^ + 0(2-u;) 



2(2 -w) 4 



(2.10) 



one can easily identify the divergent part of the polarization operator (sum of the two 
expressions ( |2.4D and (|2.5| )), with the expansion (2.6) - ( |2.9D of the divergent part of the 
effective action 

f W(diagr) = -— — i r / d^ 1 / 2 + - ~) V + 

dM,v y; 2(4vr) 2 (o; - 2) 7 y \ 60 2 V 5 6/ 



(2.11) 



that is almost the same as Eq. ( |1.7| ). The main difference is that the expression ( 2.11| ) 
does not contain surface terms jd A xg l l 2 E and Jd 4: xg l / 2 V 2 R which do not contribute 
to the propagator. 

After subtracting divergences (together with the ln[4-7r/i 2 /m 2 ] terms) and neglecting 
the 0(w — 2) terms, the polarization operator becomes 



W v ' al3 {p) 



i r i 



2(4vr) 2 1 2 



2 4^ , , P 2 ™ 2 , 3m 4 



+ 



p 4 (Ci + 5C 2 + C 3 ) - p 2 m 2 (3d + 2C 2 ) + m 4 ^ 



4 9p 4 p 2 m 2 3m 



40 



+ 



1 



2 2 

(pW 3 + pW") [p 2 (C 2 + 2C 3 ) - 2m 2 C 2 + ^ + + 2C 3 pVpV+ 
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i 2 2 

(^p-/ + ^ p y + rf a p^ + r]^p^ P a ) \p 2 (C 2 + 4C 3 ) - m 2 (d + AC 2 )-^--^- 
4 L 60 6 



* (pV-.pV) 

(pV - p 2 ^) 



Ci/ - 2dm 2 + 2p 2 C 2 - y )rf ^ + 4pVC 2 

2 

Cip 2 - 2dm 2 + 2p 2 C 2 - ^- JrT + 4pVC 2 



+ 



+ 



£m 2 



p i p ' y + p Tj rj ' — (p^p rj H + p p^rf )+ 



where 



Ci 



1/2 r ,,2 

dy In 



C 2 



1/2 /l ,, 
dy[-r-y In 



o 



(2.12) 



1/2 



are functions of the ratio p 2 /m 2 , which can be easily calculated 

C,=A, C 2 



2 r S(i-""l + l 



,4 ,4 1 

+ 



36a 2 4 36 



1 



41 



3 5a 4 6a 2 16 60a 2 3600 ' 



where 



A = l-iln 

2 - a 



4p 2 



p 2 + 4m 2 



(2.13) 



(2.14) 



For the terms without £ we can compare the expression ( |2,12| ) and the corresponding 
covariant expression 



SL = J d 4 x g 1 ' 2 { -9 G R -6 A + B X C 2 + 9 2 R 2 } 



(2.15) 



Using the bilinear expansions (2J3), ( |2.7| ), fl2.8| ), (2.9), we arrive at the following relations 
for #i, 6> 2 , #g, #A hi the momentum space (the expression ( [2.15 ) must be understood such 
that all #'s act as operators in the bilinear expansions of the corresponding covariant terms): 



(hp^ e_GP^_0_A = i_ 

2 4 4 4(4vr) 2 



C^p 4 n 9 Cim 4 p 4 p 2 m 2 3m 4 

p m C 2 H 1 1 1 

2 F 2 80 8 8 
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3 F 2 4 



1 



4(4vr) 



(Ci + 5C 2 + C 3 )p 4 p 2 m 2 {3C 1 + 2C 2 ) m 4 Ci 9p 4 p 2 m 2 3m 



+ 



80 



^ Of) \ „ 2 + ° — 1 
3~ ~ V P + T - "4(4^ 



2 2 



(2C 3 + C 2 y-2m^ 2 + f- + — 

3D o 



V _ #G _ ]_ 

2 4 ~ 4(4vr) 2 



2C 3 + ^V-(2C 2 + ^) m 2 -^-^ 
2 / V 2 / 120 12 



26 



- + 28 2 



3 ' (4vr) 2 ' 

Furthermore, it is easy to check that the ^-dependent terms can be cast into the form 



(2.16) 



(47T) 



(jfp v - p 2 rf u ) (p a p f3 - p 2 J] al3 ) 



4A 1 

2£A H =• -A + - 

. 3a 2 9 



+ 



+ 



m 



2 r 



4 L 



p r ' ^ —p rf i] ' + rj H + p prrf )+ 



1 



(2.17) 



This can be easily identified with Eq. ( |2.8[) and (p.7|), therefore the ^-dependent terms are 
associated with the J* i? 2 and J m 2 i?-terms in the effective action and contribute to the #2 
and 6c coefficients only. 

Despite the number of the equations in ( |2.16j ) exceeds the number of the unknowns, 
all these equations can be indeed satisfied. Taking into account the ^-dependent terms in 
62, we obtain 

3m 4 



; (4vr) 2 ' 



G 



m 



2(4vr 
1 



,2 1^ 



1- 

67' 



8A 2 \ 
+ 45^J 



,4 



4(4vr) 



4(4tt) 2 \15a 4 
1\2 1/^ lx/ 4 \ 2 1 1 

1 



216(4vr) 2 a 



2 „2 • 



(2.18) 



The coefficients ( 2.18] ) contain all information about the polarization operator. It will 
be checked, using the covariant Avramidi and Barvinsky - Vilkovisky technique in the next 
section. After that we shall consider the renormalization group for the coefficients of the 
vacuum action ([Q|). 
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3. The covariant derivation of the effective action 



The calculation of the polarization operator, presented in section 2, is rather complicated, 
so it is useful to verify it using a covariant method and compare the two results. To this end, 
we need to obtain the covariant effective action up to the second order in curvature. The 
one-loop effective action can be presented as the functional trace of the proper-time integral 
of the heat kernel K(s). This representation of the effective action has been obtained by 



Avramidi using summation of the Schwinger-DeWitt series [14] and by Barvinsky and 
Vilkovisky using generalized Schwinger-DeWitt technique []15|], 

In the case of a massive scalar field, the formulas of ]15| must be slightly modified 
and the one-loop Euclidean effective action (up to the second order in curvature) reads 
(compare to the formulas (1.8) and (2.1) of fj~5f| ) 

r« = -I Trlnf - V 2 i + m 2 -F + -M) = -- f°° — tr K(s) , (3.1) 
2 V 6 / 2 J s 



where P = — (£ — l/6)i? and 

,2\2-w 



where 



(..2\2-w [■ of r 

tr K{s) = K -^-^ J d t xg L ^e- am tr { 1 + sP + s 2 [R^^-sV^R^ + 

+Rf 2 (-sV 2 )R + Pf 3 (-sV 2 )R + Pf A (-sV 2 )P]} + 0(K 3 ), (3.2) 



/(r)-l + r/6 f(r) /(r)-l /( r )-l + r/6 
/l(T) = ^ ' h{T) = W + ^47 8^ ' 

f(r) /(r)-l /(r) 
/ 3 (r) = — + , / 4 (r) = — 

/(r)= [ dae^ 1 -^, r = -sV 2 . 

The integral over the proper time s is divergent and must be regularized. Below we adopt 
the dimensional regularization in the form suggested by Brown and Cassidy []l6| (see also 



and 



|17] for useful technical details). Also, inserting the exp(— sm ) factor into ( |3.2[) enables 



one to study the UV limit r/sm 2 3> 1 and the IR limit r/sm 2 <C 1 instead of the limits 



r>l and r < 1 which were investigated in [15]. The interface between UV and IR limits 
is our main subject of interest here. 

It proves useful to change the variable sm 2 = t and also denote 

u=---— 
t m 2 

After simple calculations we arrive at the following representation for the effective action 
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3 5 

+ J2 l *i- R ^ M i RIXV + J2 l r RM i R ) ' 



i=l 



(3.3) 



where 



l* x = 1 



7* — - 
h ~ 6' 



288 12 V 6/ 2 V 6 



and 
Mi 



/(fa) 

u 2 t w+l 



k = -— + -(£ 
16 2 1 



M 2 



/(fa) 



M 3 



/(fa) 



6/ ' 



Ml 



M 5 



1 



(3.4) 



The calculation of the integrals in ( |3.3j ) is quite tedious. Let us present just the result, 
using notations 



da In 



1 + a(l — a)u 



a l-a/2' 



4V 2 



V 2 - 4m 2 



(3.5) 



(indeed, this is equivalent to the previous definition Eq. ( 2.14D ), relations ( |2.10[ ) and the 
expansion (we omit, systematically, those terms which contribute only to 0(w — 2)-terms) 



m 



47T/i z 



2 \ w ~ 2 



1 + (2 - w) In 



+ .... 



The effective action can be cast into the form 



2(4vr) 



d 4 ^V 2 ^ 



-^ + ln(i^ 

2 — w V m z 



3 

+ 2 



+ 



where 



1 



60 (2 -to) 60 



^ + m(^ 

2 — w V m z 



+ 1 



+ 



+ R 



2 V s 6. 



kw(a) 



1 



2-w 
84 



+ In 



/47T/J 5 

V m 2 



+ %/(a 



Qfj.ua/3 _|_ 



i? 



2 1 

15 a 4 + 45 a 2 + 150 



(3.6) 



Ma) = 



1\ 2 A 



1 



2A / 1 



+ 



A 



A 



+ 



A 



+ 



6/ 3a 2 V s 6/ 9a 4 18a 2 144 
1 7 1/ 1\ 2 



108 a 2 2160 18 



6 



(3.7) 
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In Eq. (HH) w e used the relation C 2 = E + 2W , that was justified in |15|], including the 
terms with non-local l/V 2 -type insertions. 

It is easy to see that the above expression for the effective action perfectly corresponds 
to the polarization operator derived in the previous section. This coincidence represents 
reliable verification of the results and enables one to obtain robust conclusions concerning 
the renormalization group in curved space-time. 

Another important check of the expression ( |3.6| ) is the following. In the massless limit 
m — > we expect to meet the standard massless result, and in particular the anomaly- 
induced effective action. Taking m->Owe arrive at the singular expression for A in (|3.5|), 
On the other hand, the consistency of the conformal theory (and in particular the existence 
of the anomaly-induced effective action |L8|]) requires £ = 1/6. After we set £ = 1/6, the 
divergent ^-dependent terms cancel, and in the R[...]R sector we meet only the local term 

fW(e = l/6,m-K)) = - 12 |dW /2 * 2 +.... (3.8) 

Taking into account the well-known relation 

2 8 



we obtain the standard expression for the V 2 i?-term in the conformal anomaly 

<T ' >= ~ Tmk? v2R + ■■■ • <3 ' 9) 

where the dots stand for the E and C 2 terms which we do not discuss here. Let us remark 



J d A xg 1/2 R 2 = 12 V 2 R, 



that this test of the calculations is rather complete, because the coefficient in (3.9) depends 
on almost all terms in the expression for the effective action. Those terms which do not 
contribute to ( jO| ) directly, are verified through the cancelation of the singular (at m — > 0) 
structures. 

4. Renormalization group equations 

In the MS scheme the /3-function of the effective charge C is defined as 

Pc(MS) = l™ i/if. (4.1) 



It is easy to see that when this procedure is applied to the expressions (|3.6| ) or ( 2.18 ), the 



/?(MS)-functions for the parameters A/G, 1/G, a± t 2 are exactly the same as the usual 
ones, obtained in the covariant approach ( |1.9| ) and QL.10 ). 



The disadvantage of the MS scheme is that one can not control the decoupling of the 
massive quantum fields in an external gravitational field, hence one has to go beyond this 
scheme. The explicit calculations of the one-loop diagrams of massive fields presented in 
the previous sections enable one to apply the physical mass-dependent scheme of renor- 
malization and raise a hope to observe the decoupling. Unfortunately, this approach is not 
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universal, because the calculations have been performed through the perturbative expan- 
sion of the metric on the flat background. Also, the covariant calculation in section 3 has 
been performed in the second order in curvature approximation and can not provide an 
information beyond that obtained in section 2. Therefore, one can not be sure to obtain a 
reliable information about the non-perturbative effects. At the same time, we can always 
use the correspondence with the standard results in the MS scheme in the UV regime as a 
criterion of correctness of the results obtained in the mass-dependent scheme. 

The derivation of the /3-functions in the mass-dependent scheme has been described, 
e.g. in [13]. Starting from the polarization operator, one has to subtract the counterterm 



at the momentum p 2 = M 2 , where M is the renormalization point. Then, the /3-function 
is defined (instead of Eq. ( |4.1| )) as 

dC 



Mathematically, this is equivalent to taking the derivative —pd/dp of the formfactors in 

' ' [waft 



the polarization operator. If we apply this procedure to the formfactor of the C 2 a-term, 



the result for the /3-function in a mass-dependent scheme is 

that is the general result for the one-loop /3-function valid at any scale. As the special cases 
we meet the UV limit p 2 3> m 2 



that agrees with the MS-scheme result ( |1.10| ) and also the IR limit p 2 <C 



-i 2 4 

qlR _ 1 P , sn(P 



2 

m : 



01 ~ 1680 (4vr) 2 ' ^ + ■ ( 4 ' 5 ) 

Similar calculations for the coefficient of the i? 2 -term give 

A = l^M-^ = --^{UAA-a'A + a 2 )^- 1 -) 2 - 



n^4 dM (4vr) 2 L 8 ' \ 6/ 



a 2 -4/ 1\ , (a 2 -4) (a 2 - 12) A f V 



a 2 -4 






5A 1 


5 


+ 8 


- 6a 2 


144 


9a 4 ' 144 


108a 2 



(4.6) 



In the UV limit p 2 S> m 2 the /3-function is (in agreement with the standard result ( 1.10] )) 



% V = -^(e-s) (4.7) 



1 / 1\ 2 n ( m2 

2(4^) 2 V e ~6/ + \J 2 
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while in the IR limit p 2 <C m 2 we obtain: 



ti R 



12 (4vr) 



s 6/ 15 V s 6/ 630 



4 + 



(4.8) 



The expressions (E3)-(Em are the main results of our work. In particular, the inter- 
polation between the (3- functions in the high-energy regime ( |4.4j ), (4.7) and the /3-functions 
in the low-energy regime (4.5), (|4.8| ) exhibits the decoupling phenomenon in the vacuum 



quantum effects of the massive scalar particle. Let us remark that the correctness of the 
calculations of the /3-functions is verified by: i) precise correspondence in the UV with the 
known expressions obtained in the MS scheme; ii) great amount of cancelations of the sin- 
gular C(m 4 /p 4 ), 0(m 2 /p 2 ) and constant C(l)-terms in the expressions for j3\^{IR) in ( fO| ) 
and ( |4.Sj ); iii) correct result for the anomalous f R 2 -tevm in the massless limit. We remark 
that, due to covariance of the effective action, the expressions for the /^-functions (4.3) and 
Q4.6| ) should apply not only to the 2-point function but also to the vertex corrections. 

At the same time, taking the derivative with respect to momenta of those terms in the 
effective action, which correspond to the cosmological constant and Einstein-Hilbert term, 
brings zero result. Therefore, at this level we do not have an explicit interpretation for the 
renormalization group equations for 1/G and A/G. Unfortunately, these /3-functions can 
not be calculated in the mass-dependent scheme through the perturbative expansion of the 
metric on flat background. 



5. Discussions and conclusions 



Making the perturbative calculations in a physical mass-dependent scheme, we arrived at 
the explicit expressions for the /3-functions for the coefficients of the higher-derivative terms 
ai and a.2- The decoupling of massive degrees of freedom proceeds according to the rule 

2 

(3{IR) ~ const • fArV p 2 <m 2 , (5.1) 
V m z ) 

which has been guessed in for the cosmological constant by analogy with a similar 
phenomena in QED. However, we have discovered this law not for the cosmological constant 
but only for the higher derivative terms. 

As we have seen above, the perturbative calculations in the mass-dependent scheme 
do not reveal the /3-functions in the case of 1/G and A/G. In these sectors the UV limit of 
the mass-dependent scheme differs from that in the MS scheme. The question is whether 
these /3-functions are really zero? Prom our point of view, this is not so. The reason for 
the difference between two renormalization schemes is that actual calculations has been 
performed on the flat background, while the MS-based derivation is completely covariant. 

After all, there remains a possibility to have the scale dependence in the observables 
different from the polarization operator on the flat background. If this dependence really 
exists, it would be a non-perturbative effect with respect to the expansion (|2.1|) , since the 
perturbative effects are controlled by general covariance. It is natural to expect that these 
non-perturbative effects really take place, for otherwise there would not be correspondence 
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between the mass-dependent and MS schemes in the high-energy regime. However, it is 
not for sure that the decoupling law for these non-perturbative effects will be the same as 
for the perturbatively accessible sectors such as the higher-derivative parameters a\^- 

It is interesting to analyze further the reasons why the " perturbative" (3\ and [3\ iq are 
not seen in the mass-dependent scheme. In the perturbative approach, as for any Quantum 
Field Theory in flat space-time, the /3-functions reflect the momentum-dependence of the 
formfactors in the polarization operator. In the covariant formalism, this dependence 
corresponds to the non-local terms in the effective action of vacuum. 

We find it instructive once again to consider the QED example. Due to gauge invari- 
ance, momentum running of effective coupling constant in QED is connected with renor- 
malization of the photon propagator (the Z% renormalization constant) and the effective 
action 

has the following non-local contribution (see, e.g. (2l]]) in dimensional regularization 

for p 2 3> 77lg, 



S = —-F flu F >Ml/ + 



77 It, 

3(4tt) 2 ^ 



V 2 \ 



(5.2) 



This action can be regarded as a non-local version of the well-known Euler-Heisenberg 
Lagrangian [p2| . 

Clearly, ln(— V 2 //i 2 ) corresponds to ln(p 2 /^ 2 ) contribution in momentum space, 
therefore, the electric charge 'runs' as the momentum changes. Similarly in gravity, as 
follows from our Eq. fl3.6|) in Sect. 3 for finite part there is the following non-local contri- 
butions to the Euclidean action: 



S 



1 



2(4vr) 



d 4 x g 



1/2 



1 



Ca„ a pk w {a)C^ + Rk R (a)R 



(5.3) 



where the functions kw and k R are given by Eq. ([?j]). At high energy limit we find the 
familiar ln(— V 2 /m 2 )-type contribution (see, e.g. [^] where the high-energy limit has been 
restored from the general covariant MS renormalization group, and also |2~T| |) 



k w 



1 

60 



In 



m 2 



and 



kn 



1 



6' 



hi 



Z! 

m 2 



(5.4) 



At low energy massive particles decouple and kw and kn can be formally expanded in 
series in V 2 /m 2 to give 

kw 



840 



Z! 

2 



111 



and 



kR ~ 12 L 



1. 



1 



v ^ 6 J 15 vs 6 ; 630 



m 2 



(5.5) 



Thus, the IR running of the coupling constants a\ and a,2 is quite similar to that of QED. 
However, for the cosmological constant and the Einstein term, the situation is different. 
The reason is simple. It is clear that unlike the C 2 and R 2 terms (these two are similar 
to the Fp V F iW term in the case of QED) it is impossible to ensure renormalization group 
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running for the cosmological constant and Einstein term by inserting a certain function of 
the operator V 2 /m 2 into the action terms because this operation produces either zero (for 
the cosmological constant) or complete derivative (for the Einstein term). 

According to [Q, [5|, ^, ||, the renormalization group scaling in curved spacetime is 
defined as scaling of the metric g^ u — » e 2t g^ u . This is the key idea of renormalization group 
in curved spacetime because momentum is not defined in a general curved spacetime and, 
hence, the familiar flat space momentum scaling p — > e~ l p cannot be used. The operator 
V 2 = g^V \iS/ u scales like V 2 — > e _2 *V 2 that corresponds to the flat spacetime scaling 
p 2 — > e~ 2t p 2 . If the running for the cosmological constant and Einstein terms could be seen 
in the framework of the perturbative expansion, in the second order in h^ v it can only be 
related with the terms 

J d*x gV 2 ^(-rfysr, J d*x g^ 2 R ( - ^) 2 R 

for the cosmological constant and with the expressions 

J d*x g^ 2 R, v (-fylir, f d±x g^ 2 R(~^)R 

for the Einstein-Hilbert term. Indeed, only these terms have an appropriate scaling, are 
of the second order in curvature, and admit the ln(— V 2 /m 2 ) insertion similar to that of 

Within the perturbative scheme, one can see the renormalization group running of 
those terms only, and not of the original cosmological and Einstein-Hilbert terms. But, 
our results show that the running for these two structures doesn't take place. Probably, 
this is related to the locality of the divergences in the effective action. Another obstacle 
for the appearance of these non-local terms arises in the IR limit, where the presence of 
V 2 in denominator would violate a nice IR limit of gravity (2^, |X2|j . Indeed, the absence of 
the /3-functions for 1/G and A/G looks as an artefact of the perturbative expansion in h^ v 
rather than the fundamental property of the renormalization group in curved space-time. 
We expect that if the calculations were performed on a non-flat background, there would 
be nontrivial /3-functions at high energy in a mass-dependent scheme. This would provide 
the correspondence with the covariant MS result (|0|) at high energy and the information 
concerning the decoupling at low energy. 

Another argument supporting the existence of the running for A and G comes from 
the derivation of the effective action in the covariant proper-time cut-off regularization 
(see Appendix B), where we have found the standard /5-functions for the cosmological and 
inverse Newton constants at the UV regime and the exponential decoupling of the massive 
fields contributions to these /3-functions in the IR regime (unfortunately, the IR limit in 
the covariant proper-time cut-off regularization is not reliable). 

In conclusion, we have performed the calculations in the mass-dependent renormal- 
ization scheme and found the explicit law for the decoupling of the massive degrees of 
freedom for the higher derivative terms in the vacuum effective action. Only more compli- 
cated calculations on the non-flat background can provide a reliable information about the 
renormalization group for the cosmological and inverse Newton constants. 
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The form of the IR decoupling rule in the higher derivative sector indicates to the 
possibility of the soft transition between the stable and unstable regimes in the anomaly- 
induced inflation ||, |23|. The detailed investigation of this application will be reported 
separately. 
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Note Added. After this paper has been submitted, we learned about the work by 
Bastianelli and Zirotti 28], which is also devoted to the derivation of quantum corrections 



for the massive non-minimal scalar. The calculations of [28] are based on the worldline 
formalism, and their final result corresponds to our intermediate expressions (|2,4|), (|2.5| ) 
before the integrals from the Appendix A are taken. Technically, this integration and 
the consequent identification with covariant terms in the gravitational action are the most 
difficult parts of our work. However, the correspondence of the intermediate formulas 
provides one more useful verification of the calculations. 



6. Appendix A. The table of useful integrals 

In section 2 we have used dimensional regularization. The following integrals were relevant: 



h 



d 2w k 1 
(2tt) 2w (k 2 + m 2 )([k - p] 2 + m 2 ) ' (4tt) ? 



m 2 ) w - 2 2 F 1 {2-w,l, 



3 -p' 



k^ 



d 2w k 

(2tt) 2w (k 2 + m 2 )([k -p} 2 + m 2 ) ~ 2 
d 2w k k^k 



^ = p"h, 



(2n) 2w (k 2 + m 2 )([k - p} 2 + m 2 ) 



2 ' 4m 2 / ' 
(Al) 

(A2) 
(A3) 



T{ ] " W) (m 2 )-' 
2(47r)™ v ; 



P \ 



2Fl V- W ^2 ] ^) 



(A4) 



h = ^^P- (m 2 )- 2 

2(4tt) w v ; 



1 „ » 5 -p z \ 



2 F 1 [2-w,2, 



2 ' 4m 2 / 



(AS) 
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/ 



d 2w k 



Wk u k a 



(2ir) 2w (k 2 + m 2 )([k - p} 2 + m 2 ) 
where the parenthesis indicate to the symmetrization (e.g. D^ a ^ = l/2[D a ^ + D^ a \), 



(A6) 



2(4tt) w v ; 



Y(l-w) 2 2 / 3 -p 2 \ 

/6 = T(4^F (m) 2F A l - W ^2'^) 



(A8) 



I 



d 2w k Wk v k a kP 

(2ir) 2w (k 2 + m 2 )([k - p] 2 + m 2 ) 



h p ii p v p a p fi + 6 p^rTlPh + 3 V 1 ^ h , (A9) 



(4tt) w y ' 



(A10) 



+ 



30 ^C 2 -" ' 3 '2 ; 4^) 



T ( 1 ~ W ) (m 2 )W -l 



( 



Fi 1-«;,1,- ; 



3 -p 2 



2 ' 4m 2 



)-5^( 



l-w,2,- ; 



5 -p 2 ' 



3 -p 2 



4(4 



Here, 



2 F 1 (a,b,c;z^j = ^ 
is the hypergeometric function of z and 



n=0 



(a)n (&)n Z n 

(c) n n! 



(a) n = a(a + l)...(a + n — 1) = 



T(a + n) 



2 ' 4m 2 - 



(All) 
(A12) 

(A13) 



is the Pochhammer symbol. 



7. Appendix B. The covariant proper-time cut-off regularization 

In order to understand better our result concerning the /3-functions, let us perform an 
alternative covariant calculation using the cut-off regularization of the proper-time integral 
in the Schwinger-De Witt approach. We shall use the same proper-time integral and the 
expression for the heat-kernel as in the section 3. In some sense, the use of the covariant 
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cut-off regularization provides a smail advantage, because one can investigate not only 
logarithmic, but also quadratic and quartic divergences. 

In the covariant cut-off regularization the effective action is defined as 



fC 1 ) 



ds 



tr K(s) 



(51) 



i/n 2 



where Q is a cut-off parameter with the dimension of mass and all other notations are the 
same as in Eq. fl3.2|) . Let us perform the analysis separately for the cosmological constant, 
Einstein-Hilbert terms and for the higher-derivative terms. 

i) For the cosmological constant term, using the same change of variables as in section 
3, we obtain 



(i) 



r 



where 



ds e -m 2 s 



1 



1 



(47TS) 



2(4tt) 5 



dtt~ A e 



3 „-t 



m 



2(4vr) 



, m 



T(a, x) 



OO 

I 



dt e~ T t 



t j.a-1 



(52) 



is the incomplete gamma function (see, e.g. |25|| , formula 8.350.2). Using the expansion 

n-l 



T(— n, x) 



-l) n 



n 



r(o,x) 



E 

m=0 



{-l) m m\ 



x 



m+l 



(54) 



where F(0,x) = — Ei(— x) and 



Ei(-x) = C + ln(x) + ]T 



-x) k 



k=l 



k-k\' 



(x > 0) , 



(B5) 



for the large cut-off limit f2 2 ^> m 2 we find 



r 



(i) 



i 



4 (4vr) 2 



+ 20 2 m 2 + m 4 lnf^) 



+ 



(B6) 



The structure of divergences is the standard one. The divergences can be canceled by the 
corresponding counterterms. Let us choose them in the form 



AS = + 



1 



4 (4vr) 5 



(ft 4 + fei/x 4 ) + 2{n 2 + k 2 ^)m 2 + m 4 In 



(57) 



If the normalization condition is chosen in a simplest and natural way k\ = hi = 0, the 
finite part of the effective action will not depend on the quartic and quadratic divergences 
(see also discussion in ||). It is easy to see that the /3a can be calculated using the 
formula ( |1.9| ) and we arrive at the usual MS-scheme renormalization group equation for 
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the cosmological constant. Indeed, the same /3-function can be obtained by taking the 
logarithmic derivative of the corresponding term in (B6) with respect to the covariant cut- 
off 17 d/dVt with negative sign. Below we shall neglect the non- logarithmic divergences and 
always define /3-functions by taking the last type of derivative. 

In order to see the decoupling of the massive field at low energy, one can try to consider 
another extreme case Q 2 < m 2 . The solution (B2) gives, after we apply to it the derivative 

n(d/dn) 



Another way of getting the same formula is to use the asymptotics of the incomplete gamma 



function (see [25], formula 8.357) for large \x 
T(a,x) = x a - l e~ x [ 



(-l) m r(l -q + m) 

m=0 



+ 0(\x\~ M ) 



x m T(l - a) 



x\ 



oo 



(B9) 



and then apply the derivative Q(d/d£l) to the expression 



m 



2(4vr) 
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m 
2 



m 



2(4vr) 



^ 6 -m 2 /n 2 



(BW) 



The expression (B8) can be understood as a very fast decrease of (3\ in the region 
below the natural scale m 2 . It is tentative to consider Eq. (BIO) as a hint supporting the 
sharp cut-off decoupling for the cosmological constant ||, in contrast to the hypothesis 
Pa(IR) ~ /i 2 /m 2 • /5a (UV) which has been suggested in |l9|]. Unfortunately, the physical 
sense of the cut-off in the gravitational setting and especially the mathematical consistency 
of the Q 2 <C m 2 limit for the heat-kernel representation are not obvious. Moreover, for 
the higher derivative sector this definition is, as it will be shown below, contradictory. 
Finally, the decoupling in the lower-derivative sector of the gravitational action can not 
be understood by using the covariant proper-time cut-off, and we present it here just to 
complete the discussion. 

ii) The calculation of the divergences and the renormalization of the Einstein-Hilbert 
term goes the same way as for the cosmological constant term. The effective action is 



,(i) 



1 



CO 

1\ f ds 



6 J J (4vrs) 2 

l/Q 2 



m 2 R / 1 
2(4V)2 V~ 6 



dt t m 2 R / 1 

t 2 2(4^) 2 V ? 6 ' 
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(511) 



For Q S> m the result is 



EH 



2(4tt) 



<»"' - nr In ( ^ 



(B12) 
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The form of the counterterms and the expression for the /3-function is obvious. 



,i\ m 2 1 



« " W^"6i' <B13) 



This result agrees with the previous one Eq. (1.10) obtained in the MS scheme. In the 
low-energy limit O 2 <C m 2 we obtain 

-~ (4^)2 e • ( BU > 

Of course, this is the same type of decoupling which we already met for the cosmological 
constant case. Unfortunately, this IR limit is inconsistent, as we shall see in brief. 

iii) For the higher-derivative sector we consider, for the sake of simplicity, only the 
)/?'"' -type term. This term defines the renormalization of the C 2 -terms, and therefore 
we can expect the correspondence with the expression ( |4.3| ) in both UV and IR limits. After 
the standard change of variables we arrive at the integral representation 

oo 1 



m 2 /n 2 

where, as before, u = — V 2 /m 2 . After some algebra we arrive at the /3-function 

11 2(4vr) 26 lv 4 J V 4 6V 2 J' [ ' 

At low energies we expand exp[a(l — a)V 2 /il 2 ] into the series in V 2 /fi 2 and arrive at the 
expression 

This expressions is in disagreement with the previous one Eq. ([4.5|), which was obtained in 
the framework of the standard methods and the physical mass-dependent scheme. Hence, 
the IR limit ft 2 m 2 of the renormalization group equations in the covariant cut-off 
approach is contradictory. In fact, this contradiction shows, once again, that the physical 
interpretation of the //-dependence in gravity is quite complicated. It is possible that this 
interpretation depends on the concrete choice of the background metric and can not be 
formulated in the general form. 

The source of the difference between the expressions Q4.5D and (B17) is obvious. The 
proper-time integral representation of the effective action is formally consistent if the in- 
tegration performs from zero to infinity (see, e.g. J27||). Indeed, the lower limit of this 
integral is divergent and its UV regularization is related to the renormalization group. In 
particular, this can be seen from the UV limit Q 2 S> m 2 which is consistent in all sectors. 
But, when taking the £l 2 <C m 2 limit, we lose the link with the effective action and hence 
the corresponding results do not have much sense. 
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